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Using a three-dimensional quasiclassical technique we explore the double ionization pathways of 
a diatomic molecule driven by an intense infrared (SOOnm) ultrashort laser pulse. For intensities 
corresponding to the tunneling regime, we find that the three main ionization mechanisms have 
distinct traces when considering the sum of the momenta parallel to the laser field as a function of 
the inter-electronic angle of escape. In addition, we find that the previously observed "finger-like" 
structure in the correlated momenta of the strongly driven He is also present for strongly driven 
diatomic molecules. We show that it is mainly due to the strong interaction, backscattering, the 
recolliding electron undergoes from the remaining core. 
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I. INTRODUCTION 

Over the last few years significant progress has been 
achieved in understanding the mechanisms underlying 
the correlated electron motion in the double ionization of 
atomic systems when driven by intense laser pulses of a 
few cycle duration and infrared frequencies. Namely, for 
laser pulse intensities lO^'^-lO^^ W/cm^ — non-sequential 
double ionization — the correlated electron dynamics is 
particularly pronounced and is understood in terms of 
the rescattering mechanism The latter is a three- 
step process where first, one electron tunnels to the con- 
tinuum, then it is accelerated and, finally, it is driven 
back by the laser field to its parent ion where it transfers 
energy and liberates the still bound electron. It is now 
accepted that two are the main mechanisms leading to 
double ionization in atomic systems: a) the rescattering 
electron upon its return to the nucleus transfers energy 
to the remaining electron with both electrons ejected si- 
multaneously to the continuum-we refer to this mecha- 
nism as EI; b) the rescattering electron is ionized upon 
its return to the nucleus while exciting the remaining 
electron which is ionized at a later time at a maximum 
of the laser pulse. The latter is the so-called recoUision- 
induced excitation with subsequent field ionization mech- 
anism (REsi) [3,S,i,[ii. 

The correlated electron dynamics in the double ion- 
ization of strongly-driven diatomic molecules is a highly 
complex problem and is far from being theoretically un- 
derstood. Responding to the two-center nuclear attrac- 
tion, the correlated electron motion in the non-sequential 
double ionization regime is much more complex than 
the atomic case with new phenomena emerging such as 
bond-softening Q and alignment dependence — observed 
in molecular experiments 0, S 0|. The rescattering 
mechanism seems to be at the heart of the correlated 
electron dynamics for the case of two-electron diatomic 
molecules as was the case for the well studied He atom. 
In addition to EI and RESI, for the case of diatomic 
molecules an additional pathway is present — not previ- 



ously identified in the double ionization of He. Namely, 
upon the return of the rescattering electron to the nu- 
cleus the two electrons form a doubly excited compound 
with both electrons ionizing at a latter time p^.[Tl]| — we 
refer to this mechanism as the DE pathway. 

Given the state of the art in computational capabili- 
ties addressing the double ionization of strongly driven 
systems with three-dimensional (3-d) first-principle tech- 
niques, namely quantum mechanical ones, is an immense 
task. Currently, 3-d quantum mechanical calculations 
from first principles (ab-initio) are available for the driven 
He atom [l2j. To cope with the highly complex task 
of tackling the double ionization of diatomic molecules 
many studies use numerical quantum approaches of re- 
duced dimensionality p^ . Others use judiciously chosen 
quantum mechanical models of reduced dimensionality, 
such as the one in ref. [ll| . Yet others use semi-analytical 
quantum approaches in the framework of the so-called 
Strong-Field Approximation [l^ , not fully accounting for 
the Coulomb singularity. 

In the current paper, we study the correlated electron 
dynamics in the double ionization of diatomic molecules 
with "frozen" nuclei in the non-sequential regime as a 
function of the intensity of the laser pulse. We do so, 
using a 3-d quasiclassical technique that we have first 
developed for conservative systems, namely, the multiple 
ionization of atomic systems such as Li, Be, by single 
photon absorption [l^ . We have very recently extended 
this technique to non-conservative systems treating the 
correlated electron d yna mics of the He atom when driven 
by strong laser fields [l6| and attosccond pulses [13 . Here 
we further build on this technique by tackling more than 
one atomic centers. The advantage of this technique is 
that it is numerically very efficient. In addition, the 
method treats the Coulomb singularity with no approx- 
imation in contrast to techniques that use "soft-core" 
potentials. Fully accounting for the Coulomb singular- 
ity with no approximations is important for describing 
accurately effects such as the striking so-called "finger- 
like" structure which was recently observed [H, [l^ . This 
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structure was attributed to the strong interaction of the 
rescattering electron with the core — backscattering. In 
addition, accounting for the Coulomb singularity will be 
very important in pump-probe set-ups where VUV or 
XUV pulses are used to probe the process. For these 
latter high frequency pulses, the excursion parameter of 
the electronic motion is smaller than the atomic dimen- 
sions making it very important to incorporate effects of 
the Coulomb potential in an exact manner. 

In section II of the current paper we first briefly present 
the 3-d quasiclassical model we use and we then study the 
correlated momenta of the doubly ionized N2 diatomic 
molecule as a function of the laser field intensity. We 
identify the three different double ionization pathways 
for diatomic homonuclear molecules and show that each 
pathway leaves a distinct trace on the double differential 
probabilities of the emitted electrons. Finally, wc show 
that the striking "figure-like" structure in the correlated 
momenta identified first for atomic systems is also 

present in the double ionization of diatomic molecules. 
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FIG. 1: The correlated electron momenta of the two escaping 
electrons parallel to the polarization axis. From left to right, 
we consider three laser pulses of 10^* Watts /cm^ and a 6- 
cycle duration, 10^'' Watts /cm^ and a 3-cycle duration, and 
1.5x10^^* Watts/cm^ and a 3-cycle duration. 



switched off with a cos^ envelope, and R is the intcrnu- 
clear distance. 



III. DOUBLE IONIZATION MECHANISMS 



II. QUASICLASSICAL MODEL 

Our 3-d quasiclassical model entails the following 
steps: We first set-up the initial phase space distribution 
of the two "active" electrons in the N2 diatomic molecule. 
For intensities of the laser field that correspond to the 
tunneling regime we assume that one electron tunnels 
through the field-lowered Coulomb potential. In anal- 
ogy to using the so-called ADK quantum mechanical rate 
formula for atoms po| one can use quantum mechani- 
cal or semiclassical tunneling rate formulas for diatomic 
molecules, such as those described in refs. [2]1 . [22I . In the 
current work we use the rate provided in ref. [22l |. The 
longitudinal momentum is zero while the transverse one 
is provided by a Gaussian distribution [2^. The above 
description is valid when the instantaneous laser field is 
smaller than a threshold value which for the N2 molecule 
with first ionization energy of Ipi = 0.5728 a.u., sec- 
ond ionization energy Ip2 = 0.9989 a.u. and internuclear 
distance of 2.079 a.u. is around 0.075 a.u. To account 
for the field strengths corresponding to the over the bar- 
rier regime the initial phase space distribution of the two 
electrons is modeled by a double electron microcanonical 
distribution 

Next, we transform to a new system of coordinates, 
the so called "regularized" coordinates [25|. This trans- 
formation is exact and explicitly eliminates the Coulomb 
singularity. This step is more challenging for molecu- 
lar systems since one has to "regularize" with respect to 
more than one atomic centers versus one atomic center 
for atoms. Finally, we use Classical Trajectory Monte 
Carlo (CTMC) for the time propagation [2^. The prop- 
agation involves the "full" four-body Hamiltonian in the 
laser field with "frozen" nuclei: H — p\/2 + p\/2 — 
l/\R/2-ri\-l/\R/2-r2\ + l/\ri-r2\ + {ri+r2)-E{t), 
with E(t) the electric field, which is a cos pulse that is 



We explore the double ionization of N2 in the tunnel- 
ing regime for a laser pulse intensity of 10^* Watts/ cm? 
and 1.5 10^* Watts /cm?, respectively. For the smaller 
intensity we consider two pulses of 6 and 3-cycle dura- 
tion. Our results are obtained for laser field polarization 
parallel to the molecular axis. We note that the initial 
state described in the previous section was shown to de- 
scribe very well double to single ratios of ionization [2^ 
and accurately capture the dependence of such ratios on 
the field polarization with respect to the molecular axis. 
Future work will address including additional quantum 
mechanical effects into the model in order to accurately 
describe the dependence of, for instance, the correlated 
momenta on the angle of the polarization axis with re- 
spect to the molecular axis. From the correlated mo- 
menta of the two escaping electrons shown in Fig. [1] we 
confirm that for increasing duration of the laser pulse the 
probability for two electrons to get ionized through the EI 
pathway decreases. Indeed, the double hump structure 
in the sum of the parallel momenta distribution is much 
less pronounced for a long laser pulse versus a short one, 
see Fig. [21 as has already been pointed out in ref [l3|. As 
for the case of the strongly driven helium, the maxima in 
the parallel momentum distribution correspond to 
with 2y^l7p the maximum velocity an electron can acquire 
from its interaction with the field. Up = i?g/(4a'^) is the 
ponderomotive energy. Given the above discussion, we 
choose laser pulses of only 3-cycle duration in order for 
the traces of the EI pathway to be more pronounced and 
not be smeared out by the remaining mechanisms. 

The question wc address first is whether the three dif- 
ferent ionization mechanisms leave distinct traces when 
considering differential probabilities. In ref. [ll|, a pla- 
nar quantum mechanical model is used with the motion 
of the center of mass of the two electrons restricted along 
the polarization axis. By monitoring in time the elec- 
tron density as a function of two coordinates, that of 
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FIG. 2: The sum of the paraUel momentum distribution of the 
two escaping electrons. From left to right, we consider three 
laser pulses of 10^* Watts /cm^ and a 6-cycle duration, 10^* 
Watts/cm^ and a 3-cycle duration, and 1.5x10^'' Watts/cm^ 
and a 3-cycle duration. 



the relative motion and that of the eenter of mass of the 
electron motion, it was shown that the ionization mech- 
anisms leave distinct traces on the electron density. Our 
goal is to identify the different ionization mechanisms 
by their distinct fingerprints on the double differential 
probability of the sum of the parallel momenta as a func- 
tion of the inter-electronic angle of escape, see Fig. [3l 
To this end, we identify the recoUision time (the time 
of minimum approach of the two electrons) through the 
maximum in the electron pair potential energy. Next, 
we ideirtify the ionization time of the two electrons 
We then select three groups of trajectories. First, those 
where both electrons ionize close to the recollision time, 
which is close to (2/3-|-n)T with T the period of the field 
and label these trajectories as EI. Then, we select the 
trajectories where the recoUiding electron ionizes close 
to the recollision time, while the remaining electron ion- 
izes at a subsequent time close to a field maximum, and 
label them as RESI. Finally the trajectories where both 
electrons are ionized much later than the recollision time 
we label as DE. The double differential probability of the 
sum of the parallel momenta as a function of the inter- 
electronic angle is indeed different for each of the thus 
selected groups of trajectories, see Fig. U) 

In particular, when the electrons ionize through the 
EI pathway, they escape with large momenta and with 
inter-electronic angles less than 90°; for the trajectories 
that the electrons ionize through the RESI pathway the 
sum of the parallel components of the momenta is al- 
most zero — consistent with our experience from strongly 
driven He — and escape with angles larger than those for 
the EI mechanism but yet mainly less than 90°. Finally, 
when the electrons escape through the DE pathway the 
sum of the momenta parallel to the field is small, as for 
RESI; but the inter-electronic angle of escape is larger 
than 90°. Given the above, it is clear from Fig.[3]that for 
an intensity of 10^^ Watts/cm^ (tunneling regime) the 
longer the laser pulse the more pronounced is the contri- 
bution of the RESI and DE ionization pathways. Further 
more, for short laser pulses and intensities in the tunnel- 
ing regime for larger intensities the momenta of the two 
electrons are overall smaller, compare 10^^ versus 1.5 x 
10^^ Watts/cm^ in Fig. 01 This is consistent, with our 
finding (not shown here) that the average recollision time 
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FIG. 3: From left to right, we consider three laser pulses 
of 10^* Watts /cm^ and a 6-cycIe duration, 10^* Watts /cm^ 
and a 3-cycle duration, and 1.5x10^^ Watts/cm^ and a 3- 
cycle duration. We plot: a) (top panel) the distribution of 
the inter-electronic angles of escape binned in 14 intervals, 
180°(/ -1)/U<e < I8O714 x I with 1=1,.. .,14; b) (bottom 
panel) the sum of the parallel momenta as a function of the 
inter-electronic angle of escape. 
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FIG. 4: Top panel corresponds to a laser pulse of 10 
Watts /cm^ and a 3-cycle duration; Bottom panel corresponds 
to a laser pulse of 10^'' Watts /cm? and a 3-cycle duration. 
From left to right we show the differential probability for the 
EI, RESI and DE mechanisms respectively. 



shifts to smaller values for the larger intensity, which cor- 
responds to the electrons having smaller momenta at the 
time of their release. Finally, comparing the sum of the 
momenta as a function of the inter-electronic angle for 
the strongly driven He, see ref. [l^, with the molecular 
case shown in Fig. [3J it is clear that the RESI is more 
pronounced for the molecular case, while the DE is not 
present for the atomic case. 

In the case of the He atom, a striking "finger-like" 
structure in the correlated momenta was very recently 
observed [H, [l^ a nd different aspects of it were dis- 
cussed in refs. [3, [11]. Next, we investigate whether 
a "finger-like" structure is also present in the molecu- 
lar case and how it depends on the intensity of the laser 
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FIG. 5: Top panel corresponds to a laser pulse of 10^'' 
Watts /cm^ and a 3-cycIe duration; Bottom panel corresponds 
to a laser pulse of 10^'' Watts /cm^ and a 3-cycle duration. 
The left figures are plotted using only the trajectories where 
pi V p2 > The right figures are the same as the left 

ones except that in addition electron 2 is backscattering with 

155° <P2,,^.fP2,,,^71PlP2l <180°. 



pulse. To this end we select those trajectories for which 
electron 2 (recolliding electron) backscatters from the nu- 
cleus, inverting the direction of its velocity. That is, 

155° < P2,aft-P2,bef/\P2,aftP2,bef\ < 180° , with P2,&e//a/t 

the momentum of electron 2 just before and after the 
recoUision time. The correlated momenta of the thus 
selected trajectories, as can be seen in Fig. [SJ indeed 
correspond to a finger-like structure for the molecular 
case. As for the atomic case, we find that this struc- 
ture extends beyond the 2^yUp maximum momentum 
limit — the maximum momentum that can be acquired 
from the field. Note first that this structure persists for 
both 10^^ and 1.5 x 10^^ Watts /cm^. In somewhat more 
details in Fig. \5\ (left panel) we show the structure for 
the correlated momenta with at least one of the two mo- 
menta having magnitude greater than 2^yUp. We note 
that the trajectories shown in Fig. [5] (right panel) are 
a subset of those in Fig. [S] (left panel) and that for the 
remaining trajectories either electron 2 or electron 1 re- 
verses its velocity but with a smaller recoil angle, that 

is, 90° < Pi^aft ■ Pi,bef/\Pi,ciftPi,bef \ < 150°. It is worth 
noting that we obtain the finger-like structure in Fig.[5]in 
the first and third quadrant for electrons escaping asymp- 
totically with an angle around 50° which is larger than 
for the He case where for the recoil trajectories the two 
electrons escape almost parallel to each other. 

The distinct traces of the ionization pathways on dou- 
ble differential probabilities and the "finger-like" struc- 
ture are present for intensities of the laser pulse in the 
tunneling regime. Considering a high intensity of 10^'^ 
Watts /cm^ (over the barrier regime) we find that the 
"finger-like" structure in the first and third quadrant dis- 
appears. For this high intensity we have trajectories ion- 
izing when the instantaneous value of the field of the 
laser pulse is below the barrier and trajectories that ion- 
ize when the strength of the field is above the barrier. For 
those in the tunneling regime, we can still define a rccoUi- 
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FIG. 6: Double differential probability for a laser pulse of 10^^ 
Watts /cm? and a 3-cycle duration; Top panel corresponds to 
the trajectories that ionize over the barrier with those on the 
left corresponding to the two electrons ionizing with a time 
difference less than half a period and those on the right ion- 
izing with a time difference larger than half a period. There 
is no maximum in the electron-electron potential energy for 
these trajectories and thus no recoUision time. Bottom panel 
corresponds to the trajectories that ionize below the barrier 
(tunneling) with those on the left corresponding to electrons 
ionizing with small time difference and close to the recoUi- 
sion time; those on the right correspond to electrons ionizing 
with larger time difference with the recolliding electron ion- 
izing close to the recoUision time. For these trajectories the 
maximum of the potential energy shifts to a maximum of the 
field and thus the recoUision time is smaller compared to the 
recoUision time at 10^* and 1.5 10^* Watts/cm^. 



sion time from the maximum of the potential energy with 
the difference that the recoUision time shifts to earlier 
times at a maximum of the field compared to the recoUi- 
sion time of the smaller intensities previously considered. 
Thus, when the electrons are released close to the rec- 
oUision time their energy of release is small in contrast 
to the previously studied intensities. This is consistent 
with Fig. [HI where in the bottom panel we see that when 
the electrons ionize with small time difference and close 
to the recoUision time they escape with smaller momenta 
compared to those the electrons have for smaller intensi- 
ties (EI pathway). It can also explain that the electrons 
ionizing with larger time difference with the recolliding 
electron ionizing close to the recoUision time escape with 
higher momenta than those for smaller intensities (RESI 
mechanism); the remaining electron is not released at a 
maximum of the field resulting in higher momentum. We 
also note that the inter-electronic angle of escape is over- 
all smaller in this case due to the large strength of the 
field which pulls the electrons along its direction. Finally, 
electrons ionizing when the instantaneous field is above 
the barrier have larger momenta when the time difference 
in their ionization time is small, see top panel in Fig. [S] 
In this latter case there is no maximum in the potential 
energy and thus we can not define a recoUision time. 

In conclusion, using a 3-d quassicalssical method we 
explore the different mechanisms for double ionization of 
a strongly driven homonuclcar diatomic molecule. Wc 
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find that the ionization mechanisms leave distinct traces 
on the sum of the momenta parallel to the field polariza- 
tion of the escaping electrons as a function of the inter- 
electronic angle of escape for intensities in the tunneling 
regime. In addition, the previously observed in atomic 
systems "finger-like" structure is present for molecular 



systems as well and disappears (in the first and the third 
quadrant) for high intensities in the over the barrier 
regime. Future studies will aim to incorporate in the cur- 
rent model interference effects from the presence of the 
two nuclei and explore the "fingerprints" of interference 
on the correlated momenta of the escaping electrons. 
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